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ABSTRACT 


V/hile  the  total  number  of  theoretically  possible  quadruped  gaits 

is  quite  Jl^rge,  only  six  gaits  have  the  property  that  they  can  be  executed 

while  keeping  at  least  three  feet  on  the  ground  at  all  times.  These  gaits, 

called  creeping  gaits,  seem  to  be  well  suited  for  low  speed  locomotion 

since,  they  permit  a  quadruped  to  remain  statically  stable  during  most  of 
/ 

a  locomotion  cycle.  A  mathematical  analysis  shows,  however,  that  for 
only  three  of  the  six  creeping  gaits  is  it  possible  to  place  the  feet  of  an 
animal  or  machine  so  that  it  is  statically  stable  at  all  times.  Furthermore, 
among  these  three,  there  exists  a  unique  optimum  gait  which  mjiximizes 
static  stability.  This  gait  corresponds  to  the  normal  quadruped  crawl 
favored  by  niost  ahiimals  for  very  low  speed  locomotion. 


INTRODUCTION 


The  discrete  aspects  of  legged  locomotion  have  drawn  the  attention 
of  a  number  of  investigators  who  have  all  been  interested  in  the  problem 
of  identifying  and  class’.fying  the  natural  and  theoretically  possible  patterns 
of  footfalls  or  gaits  used  by  terrestrial  animals  for  locomotion  [  1,2, 3,4], 
More  recently,  improvements  in  technology  have  lead  to  the  development 
of  artificial  legged  locomotion  systems  and  have  raised  certain  questions 
relating  to  the  control  and  stability  of  such  machines  [  5,6,7,81.  This 
paper  is  directed  toward  the  formulation  and  solution  of  a  particular 
stability  problem  which  arises  in  connection  with  both  natural  and  artificial 
quadrupeds.  Specifically,  the  question  to  be  resolved  is  the  following: 

Among  all  of  the  theoretically  possible  quadruped  gaits,  which  ones  have 
the  property  that  the  feet  in  contact  with  the  ground  can  be  placed  so  that 
the  quadruped  is  statically  stable  at  all  times?  These  gaits  are  of  particular 
importance  for  both  animals  and  machines  since  they  represent  modes  of 
locomotion  which  can  be  employed  for  arbitrarily  slow  motion. 

The  analysis  to  follow  makes  use  of  an  extension  of  the  finite  state 
model  introduced  in  [  3  ]  to  include  the  basic  kinematic  parameters  assoc¬ 
iated  with  the  geometrical  aspects  of  a  quadruped  machine  or  animal.  The 
principal  result  obtained  is  that  only  three  quadruped  gaits  are  suitable 
for  very  low  speed  locomotion  and  that  one  of  these  possesses  stability 
properties  superior  to  the  other  two.  This  gait  is  the  typical  quadruped 
crawl  recognized  by  Muybridge  [21  as  being  the  normal  gait  employed  by 
both  natural  quadrupeds  and  human  beings  for  very  low  speed  four-legged 


locomotion. 


SUPPORT  PATTERNS  AND  STATIC  STABILITY 

Due  to  the  high  degree  of  complexity  of  natural  and  artificial  loco¬ 
motion  systems  based  upon  articulated  limbs,  it  is  necessary  to  make 
certain  simplifying  assumptions  in  order  to  obtain  some  concrete  results 
relating  to  stability.  For  the  purposes  of  this  paper,  only  steady  state 
constant  speed  locomotion  in  a  straight  line  over  a  horizontal  plane  support¬ 
ing  surface  with  the  legs  of  the  system  cycling  periodically  in  both  space 
and  time  will  be  considered.  The  legs  of  the  systems  to  be  investigated 
will  be  idealized  in  the  respect  that  the  foot  of  each  supporting  leg  will  be 
assumed  to  touch  the  supporting  surface  at  a  single  point  called  the  contact 
point.  In  addition,  it  will  be  required  that  the  force  exerted  by  any  leg  at  its 
contact  point  be  directed  into  the  t  porting  surface  and  that  no  moments 
be  applied  by  a  leg  to  the  support  plane;  i.  e. ,  locomotion  with  grasping 
feet  will  not  be  treated.  These  assumptions  and  constraints  are  formalized 
in  the  following  definition  which  imposes  some  additional  requirements. 

All  of  the  results  to  follow  are  obtained  with  respect  to  this  mathematical 
model. 


Definition  1 :  An  ideal  legged  locomotion  machine  is  a  rigid  body 
to  which  are  attached  a  specified  n'lmber,  n,  of  massless  legs.  The  length 
of  each  leg  is  arbitrarily  controllable.  Each  leg  contacts  the  supporting 
surface  at  a  point  and  can  exert  an  arbitrary  force  directed  into  this  surface. 
Arbitrary  moments  can  be  applied  to  the  body  by  any  leg  subject  only  to  the 
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constraint  that  no  moment  be  applied  to  the  supporting  surface  at  any  leg 
contact  point. 

Given  the  contact  point  for  each  supporting  leg  of  an  ideal  legged 
locomotion  machine,  a  support  pattern  T  l3for  any  phase  of  a  gait  is  defined 
as  follows : 

Definition  2:  The  support  pattern  associated  with  any  phase  of  a 
given  gait  of  an  ideal  legged  locomotion  machine  is  the  minimum  area  convex 
point  set  in  the  support  plane  such  that  all  of  the  leg  contact  points  are 
contained. 

From  this  definition  it  is  clear  that  when  m  feet  of  a  machine  are  in 
contact  with  its  support  plane,  the  support  pattern  is  a  polygon  of  not  more 
than  m  sides.  Figure  1  displays  a  sequence  of  support  patterns  for  the 
previously  mentioned  crawl  gait  TZ,  il.  The  arrows  on  this  figure  show  the 
total  motion  of  the  center  of  gravity  of  the  associated  machine  during  each 
phase  of  its  motion.  This  gait  alternates  three  and  four  feet  on  the  ground 
so  its  support  patterns  are  alternately  tri-ingies  and  quadalaterals. 

Support  patterns  are  related  to  gait  stability  for  idea.l  legged  loco¬ 
motion  machines  by  the  following  definition  and  theorem; 

Definition  3:  An  ideal  legged  locomotion  machine  is  statically  stabl ■» 
at  time  t  if  all  legs  in  contact  with  the  support  plane  at  the  given  time  remain 
in  contact  with  that  plane  when  all  legs  of  the  machine  are  fixed  at  their 
locations  at  time  t  and  the  translational  and  rotational  velocities  of  the 
resulting  rigid  body  are  simultaneously  reduced  to  zero. 

-4- 


I 


Theorem  1:  An  ideal  legged  locomotion  machine  supported  by  a 


stationary  horizontal  plane  surface  is  statically  stable  at  time  t  if  and 
only  if  the  vertical  projectitjn  of  the  center  of  gravity  of  the  machine  onto 
the  supporting  surface  lies  within  its  support  pattern  at  the  given  time. 

Proof:  If  a  machine  is  not  statically  stable  at  time  t  then,  due  to 
the  rigid  body  assumption  and  the  character  of  the  legs  of  an  ideal  legged 
locomotion  machine,  there  must  exist  an  overturning  moment  resulting 
solely  from  gravitational  acceleration.  Elementary  mechanical  consider- 
aticms  show  that  such  a  moment  exists  if  and  only  if  it  is  possible  to  draw 
a  line  in  the  support  plane  which  passes  through  at  least  one  leg  contact 
point  and  which  separates  the  vertical  projection  of  the  center  of  gravity 
onto  that  plane  from  the  contact  points  of  all  other  legs.  But  the  set  of  all 
such  lines  is  juat  the  set  of  all  tangents  to  the  support  patterti  at  time  t. 
Consequently, since  all  contact  points  are  interior  to  the  support  pattern 
and  this  pattern  is  a  convex  point  set,  a  machine  is  statically  stable  at 
time  t  if  and  only  if  the  vertical  projection  of  the  center  of  gravity  is  also 
interior  to  the  support  pattern. 

CREEPING  GAITS  AND  STABILITY  MARGIN 

While  Theorem  1  provides  a  necessary  and  sufficient  condition  for 
the  static  stability  of  an  ideal  machine,  a  real  machine  will  be  subjected  to 
various  types  of  disturbances  and  consequently  may  not  remain  upright 
even  though  the  conditions  of  this  theorem  ere  satisfied  at  all  times.  It  is 
therefore  desirable  to  introduce  a  notion  of  relative  stability  which  provides 
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s<»ne  indication  of  the  ability  of  a  machine  to  resist  disturbing  influences. 

The  following  definition  provides  such  a  measure: 

Definition  4:  The  magnitude  of  the  static  stability  margin  at  time  t 
for  an  arbitrary  support  pattern  is  equal  to  the  shortest  distance  from  the 
vertical  projection  of  the  center  of  gravity  to  any  point  on  the  boundary  of 
the  support  pattern.  If  the  pattern  is  statically  stable,  the  stability  margin 
is  positive.  Otherwise,  it  is  negative. 

Thif-  definition  leads  to  the  following  theorem: 

Theorem  2:  i'or  an  ideal  legged  locomotion  machine,  the  minimum 
number  of  legs  required  to  achieve  a  gait  with  a  strictly  positive  static 
stability  margin  at  all  times  is  equal  to  three  If  the  time  required  to 
transfer  a  leg  contact  point  to  a  new  position  is  greater  than  zero,  the 
number  is  four. 

Proof:  If  the  stability  margin  is  greater  than  zero,  then  the  area 
of  the  support  pattern  must  also  be  greater  than  zero.  This  is  possible 
only  if  the  support  pattern  has  at  least  three  sides,  implying  at  least  three 
feet  on  the  ground.  If  no  time  is  required  for  movement  of  a  leg  contact 
point  from  one  position  to  another,  then  successive  overlapping  triangles  can 
be  utilized  to  contain  a  moving  center  of  gravity  at  all  times.  However,  if 
the  leg  transfer  time  is  greater  than  zero,  v^henever  any  leg  is  lifted  from 
the  supporting  surface,  three  other  legs  must  be  in  contact  with  that  surface 
implying  at  least  four  legs  altogether.  The  sufficiency  of  four  legs  is  esta¬ 
blished  by  the  support  patterns  for  a  crawl  shown  on  Figure  1. 
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For  quadrupeds,  it  is  evident  that  all  gaits  which  satisfy  Theorem 
2  are  of  the  following  type: 

Definition  5  {Tomovic  [  9] ):  A  gait  of  an  n-legged  machine  is  a 
creeping  gait  if  every  support  pattern  involves  at  least  n-1  contact  points. 

Creeping  gaits  may  be  either  singular  or  non- singular  depending 
upon  whether  or  not  simultaneous  lifting  and  placing  of  two  or  more  feet 
ever  occurs  [  3]  .  Since  any  singular  gait  can  be  obtained  as  the  limit  of 
Come  non-singular  gait  [3],  the  analysis  to  follow  will  concentrate  on 
non-singular  gaits. 

While  it  can  be  shown  that  there  exist  a  total  of  5040  theoretically 
possible  non- singular  quadruped  gaits  [  3]  ,  only  a  small  number  of  these 
qxxalify  as  creeping  gaits.  The  next  theorem  establishes  the  number  of 
such  gaits  and  identifies  them: 

Theorem  3:  There  are  exactly  six  non- singular  quadruped  creeping 

gaits. 

Proof:  A  non-singular  quadruped  gait  always  involves  exactly  eight 
phases  since  no  two  feet  are  ever  lifted  or  placed  simultaneously  [3]. 
Since  a  creeping  gait  never  has  more  than  one  foot  off  of  the  supporting 
surface,  if  it  is  non-singular,  its  support  patterns  must  involve  alternately 
three  and  four  feet  on  the  ground.  Thus,  if  the  feet  are  nvimbered  1,2, 3,4 
as  shown  on  Figure  1,  the  succession  of  supporting  feet  is  uniquely  deter¬ 
mined  by  a  permutation  of  these  numbers  which  specifies  the  order  in 
which  feet  are  placed;  e.g.  the  crawl  of  Figure  1  corresponds  to  the 
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permutation  1432.  Since  the  first  number  in  such  a  permutation  is 
arbitrary,  it  can  always  be  chosen  to  be  equal  to  1.  With  this  choice, 
the  five  other  possible  permutations  are  1423,  1342,  1324,  1234,  1243. 

Each  of  these  is  illustrated  on  Figure  2. 

A  KINEMATIC  GAIT  MODEL 

In  [  3]  ,  a  2n-I  parameter  model  for  any  gait  of  a  n-legged  loco¬ 
motion  machine  is  introduced.  This  model,  called  a  "gait  formula", 
completely  describes  the  sequential  characteristics  of  a  gait,  but  omits 
all  of  its  spatial  properties.  The  following  definitions  extend  the  idea  of 
a  gait  formula  to  include  certain  kinematic  aspects  of  locomotion. 

Definition  6  (Hildebrand  [  11 );  The  stride  length  of  a  gait  is  the 
distance,  X  ,  by  which  the  body  of  a  locomotion  machine  is  translated 
during  any  complete  leg  cycle. 

Definition  7;  The  dimensionless  foot  position,  {x^,y.),  for  leg  i  of 
a  legged  locomotion  machine  is  a  pair  of  coordinate  values  which  specifies 
the  position  of  the  contact  point  of  any  supporting  leg.  The  origin  of  the  xy 
coordinate  system  is  the  center  of  gravity  of  the  machine.  The  x  coordinate 
axis  is  aligned  with  the  direction  of  motion  with  positive  x  directed  forward. 
The  y  coordinate  axis  is  normal  to  the  x  axis  and  is  oriented  so  that  it  is 
positive  on  the  left  side  of  the  machine.  The  scale  of  the  x  and  y  coordinate 
axis  is  chosen  so  that  X  =  1. 

Definitions:  The  dimensionless  initial  foot  position,  [y.,b.),  is  the 


value  for  the  pair  which  exists  at  the  time  leg  i  first  contacts  the 


supporting  surface  during  any  locanotion  cycle. 


The  "gait  formula"  defined  in  [  3  1  involves  two  types  of  parameters. 
The  first  of  these  is  the  duty  factor,  0^,  which  expresses  the  fraction  of 
time  that  leg  i  is  in  contact  v/ith  the  supporting  surface  during  one  complete 
cycle  of  locomotion.  The  second  is  the  relative  phase,  0^,  which  indicates 
the  amount  by  which  the  motion  of  leg  i  lags  the  motion  of  leg  1  expressed 
as  a  fraction  of  the  time  required  to  complete  one  cycle  of  locomotion. 
These  parameters  together  with  the  initial  foot  position  variables  are 
combined  in  the  following  definition: 

Definition  9:  A  kinematic  gait  formula,  k,  for  an  n-legged  loco¬ 
motion  machine  is  the  (4n-l)-tuple 


k  = 


(^l*  02*  •  *  •  1*^2’  * 


V  ^’^2’ 


n 


..0  ) 

n 


(1) 


The  following  theorem  establishes  the  extent  to  which  a  kinematic 
gait  formula  determines  ihe  motion  of  a  locomotion  machine. 

Theorem  4:  A  kinematic  gait  formula  completely  specifies  the 
position  of  the  supporting  feet  and  thevertical  projection  of  the  center  of 
gravity  of  an  ideal  legged  locomotion  machine  in  both  time  and  space  up 
to  a  multiplicative  factor  of  t  in  time  and  X  in  space. 

Proof:  If  t'  stands  for  a  time  variable  measured  in  arbitrary  units 
and  T  is  the  cycle  time  [3]  required  for  one  complete  locomotion  cycle 
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measured  in  the  same  units,  then  a  dimensionless  time  variable  can  be 
defined  as  t  =  t'/r.  If  t  =  0  is  chosen  to  correspond  to  the  time  when  leg  1 
first  touches  the  supporting  surface,  then  since  is  a  spatially  normalized 
variable 

Xj  (t)  =  Yj  “  t  5  ® 

Since  leg  1  does  not  touch  the  ground  during  B  <  t  <  1,  (t)  is  undefined 

for  this  period  of  time.  Since  t  is  normalized,  the  period  of  x^(t)  is  equal 
to  1  so  (2)  determines  Xj(t)  for  all  t  such  that  leg  1  is  in  a  supporting  phase. 
For  any  other  leg  i, 

x.(t}  =  Y-“  0-  <  t<  0.  +  6.  (3) 

I  11  1  ~  1  1 

and  x,(t)  is  undefined  for  0,+  8.  <  t  <  0.+  1.  Again,  since  the  period  of  leg  i 
is  also  equal  to  1,  (3)  determines  the  dimensionless  leg  oosition  for  any 
supporting  leg. 

Since  the  motion  of  the  machine  is  entirely  along  the  x-axis  and  the 
leg  positions  are  required  to  be  periodic  in  space,  y.=  6.  whenever  x.  is 
defined.  Thus,k  determines  (x,,y.)  for  all  supporting  legs  for  all  t.  More¬ 
over,  if  is  dimensionless  location  of  the  vertical  projection  of 

the  center  of  gravity  of  the  machine  relative  to  its  initial  positon  on  the 
support  plane,  then 

so  the  center  of  gravity  location  is  also  determined.  Finally,  multiplication 
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of  an/  dimensionless  pair  W  ^  produces  positions  in  the  units  of  X, 


OPTIMUM  KINEMATIC  GAIT  FORMULAS 


The  crawl  gait  shown  on  Figure  1  is  described  by  the  kinematic  gait 


formula 


.  -  .  .  n  .  1  11  113  1 

V  12’  12’  12  ’  12  ’  '  ’  4’  ‘4’  4’  ”4*  2’  4’  4  ^ 


While  examination  of  this  figuie  shows  that  this  gait  possesses  a  positive 
stability  margin  at  all  times,  it  is  possible  that  some  other  kinematic  gait 
formula  would  produce  a  more  stable  crawl.  It  is  useful,  therefore,  to 
attempt  to  find  optimtim  gait  formulas  for  all  six  of  the  creeping  gaits  of 
Figure  2.  The  following  definitions  make  this  notion  precise. 

Definition  10:  The  dimensionless  longitudinal  stability  margin  at 
time  t  for  an  arbitrary  support  pattern  associated  with  a  kinematic  gait 
formula,  k,  is  denoted  by  s(t,k)  and  is  equal  in  magnitude  to  the  shortest 
distance  along  the  x  coordinate  axis  from  the  center  of  gravity  of  a  loco¬ 
motion  machine  to  an  edge  of  its  support  pattern.  If  the  pattern  is  statically 
stable,  the  stability  margin  is  positive.  Otherwise,  it  is  negative. 

Definition  11:  A  kinematic  gait  formula,  k  ,  for  a  given  gait 
matrix  [3],  G,  is  longitudinally  optimum  if 


max  min  8(t,k)  =  min  8(t,k*)  =  s*(G)  (6) 

keK^te[0,l]  t€[0,l] 

where  K  is  the  set  of  all  gait  formulas  which  imply  G  [3].  The  function 

(j 


The  connection  between  s*(G)  and  static  stability  for  quadruped 
creeping  gaits  is  established  by  the  following  theorem: 

Theorem  5:  Let  G  be  the  gait  matrix  for  a  quadruped  creeping 
gait  and  let  be  the  set  of  all  kinematic  gait  formulas  which  imply  G. 

Then  there  exists  a  k  e  K _  such  that  the  machine  associated  with  k  possesses 
a  strictly  positive  static  stability  margin  at  all  times  if  and  only  if  s*{G)  >0. 
Moreover,  when  such  a  k  exists,  it  satisfies  the  constraints:  6^  >0,  6^  >  0, 

6-  <0,  6  .  <  0. 

Z  4 

Proof:  Suppose  the  condition  on  the  6  variables  is  satisifed.  That 
is,  the  m.achine  has  two  feet  on  each  side  of  its  centerline.  In  this  case, 
the  projection  of  its  center  of  gravity  is  in  the  interior  of  its  support  pattern 
and  it  has  a  strictly  positive  static  stability  margin  if  and  only  if  its  longi¬ 
tudinal  stability  margin  is  strictly  positive.  But  this  is  satisfied  at  all 
times  if  and  only  if  s*(G)  >0.  Now  suppose  the  condition  on  the  6  variables 
is  not  satisfied.  Then  there  must  exist  sane  phase  of  the  gait  such  that  the 
center  of  gravity  of  the  machine  lies  either  on  the  right  or  left  boundary  of 
the  support  pattern  or  entirely  outside  of  it.  Consequently,  independently 
of  s’f'(G),  such  a  gait  will  fail  to  have  a  strictly  positive  static  stability 
margin  at  some  time  in  any  complete  locomotion  cycle.  Thus,  k  exists  if 
and  only  if  8*(G)  >  0  and  6^  >0,  6^  >  0,  6^  <0,  6^  <  0. 
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OPTIMIZATION  OF  THE  CRAWL  GAIT 

Examinat'.on  of  Figure  1  shows  that,  due  to  the  forward  motion  of 
the  center  of  gravity  of  a  locomotion  machine,  the  minimum  value  of  the 
longitudinal  stability  margin  associated  with  any  support  pattern  of  any 
gait  occurs  either  at  the  moment  the  pattern  is  established  or  at  the  end 
of  the  interval  in  which  the  pattern  exists.  In  addition,  if  the  minimum 
occurs  at  the  beginning  of  the  interval,  it  is  associated  with  a  line  at  the 
rear  of  the  support  pattern  while  if  it  occurs  at  the  end,  it  is  associated 
with  a  line  at  the  front  of  the  pattern.  Thus,  since  a  total  of  eight  support 
patterns  comprise  a  non-singular  crawl  gait,  there  are  eight  critical  times 
at  which  the  minimum  longitudinal  stability  margin  may  occur.  These  are 
all  indicated  in  the  following  set  of  eight  linear  equations  which  establish 
the  position  of  the  critical  feet  at  the  times  indicated. 


Leg  4  lift-off 

*3  =  ^3+ *3-  ®4 

(7) 

■■z  “  ''2*  *2-  *4-  «4 

(8) 

Leg  4  touch-down 

*1  =  '<1- 

(9) 

*2  '  ''2^  »2-  **4-  ‘ 

(10) 

Leg  2  lift-off 

*3  “  ''3*  *3"  *2‘  ®2 

(11) 

*4  '  ^'4^  *4-  "z-  ®2*  ‘ 

(12) 
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In  order  to  calculate  the  longitudinal  stability  margin  at  each  of  the 
above  times,  it  is  necessary  to  assume  some  values  for  the  y-coordinate 
of  each  leg.  Since  all  natural  quadrupeds  possess  right-left  symm.etry,  it 
will  be  c''*'venient  to  make  thi«  assumption  here.  In  addition,  to  simplify 
the  analysis  still  further,  it  will  also  be  assumed  that  the  spacing  between 
the  front  legs  in  the  y  direction  is  the  same  as  for  the  rear  legs.  Thus, 
the  class  of  gait  formulas,  K_,  to  be  considered  satisfy  the  constraint 

Li 
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! 


5i  =  6^  =  -5^  =  .6^  =  6  >  0 


(23) 


With  this  constraint, the  longitudinal  stability  margin  is  independent  of  6. 
In  particular,  for  the  crawl  gait,  from  Figure  1  and  eq.  (7)  and  (8),  at 
t  =  tj,  the  horizontal  distance  from  the  center  of  gravity  to  the  rear  of  the 
support  patte’-n  is  just 


£,W  =  =  i  {-V2-  Y3-  03-  20^+  ) 


At  t  =  t^,  the  distance  to  the  front  of  the  support  pattern  is 


X  X 

£,W  =  -L-1  =  i  ^  02-  20^-  1  1 


(24) 


(2  5) 


Consequently 


min  x(t,k) 


=  min  f.  (k) 

ieEl.Zl  ^ 


(26) 


Continuing  this  analysis  for  the  other  support  patterns,  and  associating 
f^(k)  with  conditions  at  t  =  t^,  from  eq.  (11)  through  (22),  the  appropriate 


functions  are; 


f3(k) 


■^3  ■''4  1 


■=  li-'s- V  «3-«4£202£«a-i  l 


^1  ^  1 

f  (k)  =  — - -  =  — -fv+Y+fi-Z®  1 

4'  '  2  2  ^4  4  2  J 


^  =  H-^1- 


(27) 


(28) 


(29) 
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(30) 


f6(k)  = 

f-{k)  = 


xi+x^ 


=  j  \"Y3"  Y^"  ^4^  ^^3  ) 


(31) 


fsik)  = 


X.,  +  x. 
<.  3 


=  |■{Y2+Y3+^2+e3-2'^ 


(32) 


Finally,  by  analogy  to  (26) 


min  8(t,k)  =  min  f.  (k)  =  f  (k)  (33) 

t€  [0,i;j  ie{l,2,...8}  ^  ° 


and  k*  is  found  by  maximizing  the  criterion  function  f  (k)  over  K„. 

————————  o  (.1 

Referring  again  to  Figure  1,  it  is  evident  that  increasing  the  fraction 
of  time,  0.,  that  leg  i  spends  in  contact  with  the  supporting  surface  never 
reduces  the  static  stability  margin  of  a  gait.  Thus,  for  any  gait,  f^(k)  is 
non-decreasing  in  6.  and  at  least  one  k*  will  involve  0^=  1  for  all  i.  But 
this  is  impossible  for  a  non-singular  gait  [3].  To  avoid  this  difficulty, 
the  further  restriction  will  be  placed  on  that 

0j  =  02  =  03  =  6^  =  0  (34) 

i.  e. ,  only  regular  gaits  [3]  will  be  considered.  The  criterion  function 

f  (k)  will  then  be  optimized  for  a  fixed  value  of  0 . 
o 

Even  after  the  parameters  0^^  are  fixed  according  to  (34),  f^(k) 
remains  a  nonlinear  function  of  seven  variables  to  be  minimized  subject 
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to  the  inequalities 


>1-3 

(3  5) 

>  0^+  (1  -  3) 

(36) 

03  >  02+ 

(3  7) 

which  are  dictated  by  the  sequence  of  support  patterns  involved  in  a  Ci*awl. 
Problems  of  this  type  fall  within  the  domain  of  nonlinear  programming  [  10] 
and  are  typically  very  difficul':  to  solve.  In  this  particular  case,  however, 
a  solution  has  been  found  by  means  of  relatively  straightforward  iterative 
procedure.  Rather  than  recording  the  details  of  this  calculation,  the 
results  are  presented  as  a  theoiem; 

Theorem  6:  If  6^=  6^-  ”^4“  1  >  0  >  3/4,  any  longi¬ 

tudinally  optimum  kinematic  gait  formula  for  a  regular  crawl  gait  must  be 
of  the  form 


k  =  3.6,9, a+|,a+|,-a+|, -a+|,6, -6. 6,-6, =  k*  (38) 

where 

a  >  I  ,  6  >  0  (39) 


Furthermore,  such  a  crawl  is  a  singular  crawl  with  the  gait  matrix  [3l 


G 


/  0  0  0  o\ 

/  0  0  0  i  \ 
j  0  1  0  0  I 

1  0  0  0  0  I 
\  0  0  1  0  / 
\  1  0  c  0/ 


(40) 
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The  longitudinal  stability  margin  for  this  gait  is  given  by 

s*(G)  =  5-3/4  (41) 

Proof:  First  of  all,  it  is  clear  that  the  constraints  imposed  by 
(3  5),  (36),  and  (37)  are  satisfied.  Then,  substituting  k*  into  (24),  (2  5), 
and  (27)  through  (32), 


tl(k*)  =  s-| 

(42) 

=  a.f4 

(43) 

=  a-|.  i 

(44) 

£^(k»)  =  e-  1 

(45) 

£5(k»)  =  8-1 

(46) 

=  a  -  1  £  i 

(47) 

fylk,,  =  a  -  I  +  i 

(48) 

£g(k<.)  =  8-1 

(49) 

Thus,  the  criterion  function  f  (k*)  is  determined  either 

o 

by  f^,  f^,  f^,  and 

fg  or  by  f^,  fg,  and  f^.  To  determine  which  of  thes 

e  cases  applies,  it 

is  necessary  to  consider  the  allowed  range  of  values  for  the  arbitrary 
constant  a.  A  lower  bound  on  a  can  be  obtained  from  the  requirements 

Xi(t)  >  x^(t)  (50) 

and 

X2(t)  >  x^{t)  (51) 
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imposed  by  the  necessity  of  avoiding  interference  between  front  and  rear 

feet.  From  Figure  1  it  can  be  seen  that  the  minimum  value  for  x,  (t)  -  (t) 

occurs  at  t  =  t-  at  which  time 
o 

x^(t)  -  x.^{t)  =  Yj-  **3  "  Y3  =  2a  -  0  >  0  (52) 

so 

a  >  I  (53) 

and  for  i  =  2 , 3 , 6,  7 

f.(k*)  >  ^  (54) 

Solution  of  (51)  yields  identical  results.  On  the  other  hand,  for  i  =  1,4,  5,8, 
since  0  <  1 , 

f.(k«)  <  i  (55) 

It  follows,  therefore,  that  the  longitudinal  stability  margin  is  determined 

by  f ,  I  f ,  I  f  _»  and  f„  for  all  allowed  values  of  a. 

14  5  8 

Since  (42),  (45),  (46)  and  (49)  all  yield  the  same  values  for  the  k 
given  by  (38),  this  value  is  optimum  unless  there  exists  some  vector,  Ak, 
such  that 

f  (k*  +  Ak)  >  f  (k*)  (  -r  ,. 

o  o 

This  in  turn  is  possible  if  and  only  if  f^,  f^,  f^,  and  fg  can  all  be  increased 
simultaneously.  To  determine  whether  or  not  this  is  possible,  from  (24), 
(28),  (29),  and  (32),  substitution  of  k*  r  Ak  for  k*  results  in  the  reladons^nps 
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:  } 
•  j 

;  i 

I  i 

I 


=  ^1  Ay^+  Ay^+  A0^-  2A02  f  (58) 

Af^  =  i  I  -Ay  Ay^-  A0^+  2A0^  '■  (59) 

Afg  =  I  {  AY2+  AY3+  A02+  A0g  I  (60) 

where 

Af.  =  f.(k=i<  +  Ak)  -  f.(k*)  (61) 


Since  0^=  8  in  (3  9),  the  left  hand  side  of  (38)  is  satisfied  as  an  equality. 
Consequenily,  it  is  necessary  in  (57)  through  (60)  that 

A0g  <  0  (62) 

Likewise,  since  0  =  and  0  =  B  -  in  (38),  (36)  is  also  satisfied  as  an 

C0  Cd  ^  Ct 

equ.^lity  with  the  consequence  that 

A0.,  =  A0.+  e  ,  €  >  0  (63) 

— 

must  also  be  satisfied  in  the  evaluation  of  (57)  through  (60). 

Substituting  (6.'  }  into  (57)  through  (60)  and  requiring  that  Af.  >  0, 
i  =  1,4,  5,8,  results  in  the  expressions 


-Ay_  -  jY,  -  A0,  -j  A0,  -  e 

C  l  3  4 

>  0 

(64) 

Ay,  ~  6v.  +  A0-  -  2c 

Ini 

>  0 

(65) 

-Ay,  “  +  2A0-  -  A0 

1  4  3  4 

>  0 

(66) 

Ay,  ’  Ay,  +  A0,  +  A0  -J  € 

^  i  i  4 

>  c 

(67) 
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Addition  of  these  four  equations  produces  the  relation 


2A0^  -  2e  >  0 


A03  >  €  >  0 


But  this  is  a  contradiction  since  (62)  requires  <  0.  Consequently,  no 
Ak  exists  which  simultaneously  increases  and  fg  so  (38)  is 

optimum.  The  G  matrix  associated  with  (38)  can  be  constructed  from  the 


algorithm  given  in  [  3  ]  . 


To  show  the  uniqueness  of  (38),  it  is  sufficient  to  equate  (64)  through 
(6V)  to  zero.  In  that  case,  (69)  becomes 


A03  =  e  ^  0 


In  view  of  (62),  this  implies 


A03  =  €  =  0 


which  in  turn  requires  that  60^  and  A02  also  equal  zero  in  order  for  (64) 

through  (67)  to  be  satisfied  as  equalities.  Thus, there  is  no  Ak  ^  0  such 

that  f  (k*  +  Ak)  =  f  (k=!=)  and  (38)  is  therefore  unique, 
o  o 


OPTIMIZATION  OF  OTHER  CREEPING  GAITS 


In  the  appendix  to  this  paper,  equations  alalogoas  to  (24)  through  (32) 
are  developed  for  each  of  the  f.ve  other  creeping  gaits.  These  equations  can 
be  subjected  to  an  optimization  process  in  the  same  way  as  the  crav-l 


was  treated. 
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The  results  of  such  an  optimization  are  detailed  in  the  appendix  and 
summarized  in  Table  I.  This  table  provides  the  basis  for  the  following 
theorem: 

Theorem  7:  For  any  specified  value  of  the  duty  factor,  B.  in  the 

3 

range  ^  <  B  <  It  the  optimum  dimensionless  longitudinal  stabili  /  margin 
for  a  regular  crawl  is  greater  than  for  any  other  regular  gait.  Moreover, 
there  are  only  two  other  non-singular  creeping  gaits,  1342  and  1243,  which 
possess  kinematic  gait  formulas  with  a  strictly  positive  value  for  their 
minimum  longitudinal  stability  margin  for  any  B  <  1. 

Proof:  Table  I  proves  the  first  part  of  the  theorem  and  the  second 
part  with  respect  to  regular  gaits.  For  irregular  gaits,  it  is  sufficient  to 
note  that  since  the  minimum  longitudinal  stability  margin  is  non-decreasing 
in  B^.  if 

p  =  max  B.  (72) 

m  1 

then  the  regular  gait  with  B  =  6  will  be  at  least  as  stable  as  the  associated 

m 

irregular  gait.  It  follows,  therefore,  that  since  the  stability  margin  of 
regular  kinematic  gait  formulas  for  the  gaits  1324,  1432,  and  1234  is 
bounded  by 

s*(G)  =  B-1  <  0  (73) 

this  bound  also  applies  to  the  irregular  gaits. 
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DISCUSSION  OF  RESULTS 


mf. 


sr 
~ 
i  ^ 


Ip 


I 


r  I 


The  information  contained  in  Table  I  is  in  dimensionless  form.  In 
absolute  form,  in  any  units,  since  s*(G)  is  normalized  to  X,  the  denormal- 
ized  optimum  stability  margin  is  given  by 


X  =  Xs*(G) 
s 


(74) 


This  shows  that  it  is  desirable  to  make  X  as  large  as  possible  consistent  with 

satisfying  the  requirement  for  no  interference  between  legs  expressed  by 

(52).  Since  the  variable  a  in  this  expression  is  normalized  to  X  while  S  is 

not,  if  X  is  the  denormalized  value  of  a,  (52)  becomes 
a 


2x  -  XS  >  0 
a  — 


(7  5) 


so  the  maximum  value  for  X  is 


2x 


m  B 


The  corresponding  value  for  a  is 


(76) 


m 


X 


m 


B 

2 


(77) 


so  from  (38),  the  unique  regular  crawl  gait  which  maximizes  the  denormal¬ 
ized  minimum  longitudinal  stability  margin  is  defined  by 


k  =  (  8, 8, 8, 8, 6, 6, 0, 0,  6,  -  6,  6,  - 6,  4 ,  B,  B  -  4 
V  2  2 


(78) 


which  involves  only  the  two  parameters  6  and  B*  This  gait  is  illustrated  by 
Figure  3. 
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Examination  of  Muybridge's  photographs  [2]  shows  that  Figure  3 
corresponds  closely  both  in  space  and  time  to  the  normal  slow  quadruped 
walk.  The  only  important  exception  is  that  since  (78)  requires  that  each 
rear  foot  be  placed  in  the  footprint  of  the  foot  ahead  of  it  at  precisely  the 
time  the  front  foot  is  lifted,  living  quadrupeds  employ  8  <  3/4  so  that  the 
front  foot  is  lifted  a  short  time  before  the  rear  foot  is  placed  in  order  to 
provide  the  necessary  clearance  betv/een  feet.  The  rear  foot  is  sometimes 
placed  somewhat  ahead  of  the  front  foot  footprint  to  compensate  for  the 
early  lifting  of  the  front  foot  [  1  ]  .  Apparently,  for  moderately  low  speed 
locomotion,  this  solution  to  the  foot  interference  problem  is  more  satis¬ 
factory  for  the  majority  of  living  qtaadrupeds  than  the  alternate  solution  of 
shortening  the  leg  stroke  of  the  optimum  crawl  defined  by  (78).  It  appears 
that  a  true  crawl  is  generally  used  only  for  grazing  and  for  other  very  low 
si^eed  activities  [  1  ]  . 

Hildebrand  [  ll]  has  observed  that,  apparently,  no  living  qxiadruped 
uses  the  gait  1324  (gait  number  5  in  his  classification  scheme)  and  has  con¬ 
jectured  that  this  might  be  due  to  the  poor  stability  properties  of  this  gait. 
The  analysis  presented  here  confirms  this  hypothesis.  It  seems  likely  that 
the  gaits  1432  and  1234  are  also  unused  for  the  same  reason.  The  authors 
are  unaware  of  any  observations  of  the  gaits  1342  and  1243  in  use  by  a  living 
quadruped  even  though  these  gaits  are  stable  for  a  sufficiently  large  duty 
factor.  It  may  well  be  that  th>  optimality  of  the  crawl  has  caused  it  to  be 
adopted  as  the  unique  gait  utilized  by  natural  quadrupeds  for  very  low  speed 
locomotion. 
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SUMMARY  AND  CONCLUSIONS 


The  results  obtained  in  this  paper  show  that  there  exists  a  unique 
solution  to  the  problem  of  choosing  a  gait  to  optimize  the  stability  of  low 
speed  quadruped  locomotion.  This  gait  corresponds  essentially  to  the  low 
speed  gait  preferred  by  most  natural  quadrupeds.  The  analysis  excludes 
certain  types  of  locomotion  due  to  assumptions  on  the  geometry  of  the 
machine  or  animal.  In  particular,  the  machine  is  not  allowed  to  drag  a 
stabilizing  tail  as  kangaroos  and  some  legged  reptiles  do  [  12  ]  ,  nor  is  it 
permitted  to  overlap  its  feet  as  would  be  possible  if  the  front  pair  of  legs 
were  spaced  more  widely  or  more  narrowly  than  the  rear  pair. 

While  the  conclusions  reached  in  this  paper  are  obtained  with  respect 
to  an  idealized  model,  the  model  appears  to  be  sufficiently  close  to  reality 
to  permit  qualitative  extrapolation  to  real  locomotion  systems.  Obviously, 
such  characteristics  of  real  systems  as  the  non-zero  mass  of  legs  could 
be  included  in  a  more  complicated  analysis.  Such  analysis  should  probably 
be  carried  out  for  particular  machines  whenever  it  is  important  that  the 
greatest  possible  degree  of  static  stability  be  achieved. 

The  type  of  analysis  presented  here  could  be  readily  extended  to 
a  study  of  insect  gaits.  It  would  be  interesting  to  determine  if  iisects  also 
utilize  optimally  stable  modes  of  locomotion.  It  appears  that  stability 
analysis  of  biped  gaits  requires  the  use  of  dynamic  models  which  ar**  neces¬ 
sarily  more  complicated  than  the  kinematic  model  defined  in  this  paper,. 

The  same  is  true  of  higher  speed  quadruped  gaits  which  typically  contc-.Li 
statically  unstable  phases  [  1]. 
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Figure  2:  Quadruped  creeping  gaitp 
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APPENDIX 


Stability  Margin  Optimization 


Gaits  1  3  4  2  and  12  4  3 


The  foot  position  equations  for  gait  1  342  are  as  follows: 


Leg  3  lift-off 

^1  “  "^1  ^ 

(79) 

"^4  =  \  +^4-'P3-^ 

(80) 

^2= 

Leg  3  touch-down 

X,  =  Vj-Pj 

(81) 

*2  = 

(82) 

Leg  4  lift-off 

(83) 

"2  =  ''2+'P2-'f4-« 

(84) 

^4= 

Leg  4  touch  down 

*1  =  ''r‘»’4 

(85) 

*2  =  ^2*'f2-'l>4-' 

(86) 

5* 

Leg  2  lift-off 

X  =  Y.+  cc,-cp,-P  +1 

(87) 

’^4  =  V'f4-V® 

(88) 

’  5’ 


Leg  I  touch-down 


Xj  =  Y]  -cp^  (89) 

=  y  4^  ^4-^1 

Ths  equations  for  t  and  t  are  the  same  for  all  gaits  and  are  given  by 
7  o 

(19),  (20),  (21),  (22).  Combining  the  foot  position  equations  appropriate!  / 
leads  to  the  s<-ability  margin  equations- 


9  = 

2 

1 

2  ' 

(91) 

‘2  = 

Xj+X^ 

2 

j_ 

2 

»  (92) 

■’^2-^3 

2 

1 

2 

(-V2-Y3-92-P3+V26-1] 

>  (93) 

1 

\ 

(\+V*2-^4-‘ 

)  (94) 
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The  other  two  stability  margin  equations  are  again  the  same  for  all  gaits 
and  are  given  by  (31)  and  (32). 

Examination  of  the  above  equations  shows  that  the  unique 
optimvim  kinematic  gait  formula  for  this  gait  is  given  by 


k=  (8,3,p,  3.  Yj.  Y^,  'V,-5  +  -  ,  -Yj-*-  33-  j,  6.  -6.  6rS,P,  33-2,  23-1) 


(97) 


For  this  gait  formula,  if  Y^  and  Y^  are  sufficiently  large,  then  the  long'- 
tudlna’  stability  margin  is  determined  by 


s*(0).  t,  =  =  |e-| 

Since  reflection  of  gait  1  342  about  an  axis  parallel  to  the  centerline  of 
'.ts  support  patterns  produces  gait  1-54  (Figure  21,  (98)  applies  equally 
well  to  either  gait. 


Gaits  1234  and  1432 


r'or  gait  1432,  the  stability  margin  equations  are: 
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(102) 

f-  = 

0 

(103) 
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!  \ 

'6  "" 

2  ^ 

(104) 

and  f  and  f„  are  again  given  by  (31)  and  (32).  An  optimum  kinematic  gait 
7  o 

formula  for  this  set  of  equations  is: 

k  =  (S,  B,  3,  3,  a+p/2,  a-f-3/2,  -ai-3  /2,  -a+3/2,  6,  -6,  6,  -6,  3/4, 

3-1/4.3-1/2) 

(105) 

witn  this  set  of  gait  parameters,if  the  variable  a  is  a  sufficiently  large 
positive  number,  then  the  stability  margin  is  determined  by 

s*(G)  =  =  3 -i  (106) 


-in- 


4 


The  k  given  by  (105)  is  not  unique  since  if  Acp^  =  ZA^  and 
^2  ~  ‘^^3*  xinchanged.  Since  gait  1234  is  a  reflection 

of  gait  1432  about  a  longitudinal  axis,  (106)  applies  to  both  gaits. 

Gait  J  324 

For  gait  1324: 
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(114) 

The  corresponding  unique  optimuna  kinematic  gait  formula  is: 

k  =  (0,5,  0,0,a+0/2,a+0/2,a+b/2,  -a+0/2,  -a+P/2,  6,-6.  6.  -6,  l/2,e-l/2.e) 

(115) 

With  this  formula,  if  the  variable  a  is  sufficiently  large  and  positive, 
s*(G)  =  f^  =f^  =  f3  =  fg  =  0-l 


(116) 
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This  paper  examines  all  of  the  theoretically  possible  quadnq>ed  gaits.  Only 
six  gaits  can  be  executed  vfaile  three  feet  are  on  the  ground  at  all  times. 
These  six,  called  creeping  gaits,  pennit  a  quadsrtqped  to  remain  statically 
stable  diffing  most  of  a  locomotion  cycle.  The  gaits  were  analyzed  by  means 
of  a  2n-1  parametOT  model  with  extensions  that  account  for  basic  kinematic 
parameters  associated  with  geometric  aspects  of  a  quadruped  machine  or 
animal.  The  model  considers  only  s'^  -'dy  state,  constant  speed  locomotion 
in  a  straight  line  over  a  horizont__  plane  supporting  surface  with  the  legs 
of  the  system  cycling  periodically  in  both  space  and  time.  Only  three 
of  the  six  gaits  enable  the  quadruped  to  be  statically  stable  at  all  times. 
Of  the  three,  one  is  a  \iniqae  optimum  gait  with  .Tiaxinium  static  stability. 

The  gait  corresponds  to  the  normal  qaadnqped  crawl  facvored  by  most  ani  mal  s 
for  very  low  speed  locomotion.  The  model,  though  idealized,  sq)pears 
sufficiently  close  to  reality  to  permit  qualitative  extr^lation  to  real 
locomotion  systems.  Stability  analysis  of  biped  gaits  requires  the  use 
of  dynamic  models  that  will  necessarily  be  more  congplex  than  the  kinematic 
model. 
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